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ON QUATERNIONS AND THEIR GENERALIZATION AND THE HISTORY 
OF THE EIGHT SQUARE THEOREM. 

By L. E. Dickson. 

1. Objects of the paper. We shall present the history of the generaliza- 
tions to four and eight squares of the familiar formula 

(1) (a 2 + 6 2 )(a 2 + /3 2 ) = r 2 + s 2 , r = aa - bp, s = a/3 + 6a, 

and an elementary exposition of Hurwitz's proof that such a formula 
holds only for 2, 4 or 8 squares. For these three cases we shall show 
that the formula admits of a simple interpretation concerning the norms 
of numbers which are ordinary complex numbers, quaternions or numbers 
of Cayley's algebra with 8 units. No knowledge of quaternions or the 
latter algebra will be presupposed, but their more fundamental algebraic 
properties will be developed in detail. 

A clear exposition will first be given (§§ 1-5) of the main results of 
our subject. This will be followed (§§ 6-28) by an account of its history, 
which is believed to omit no paper on the eight square theorem and its 
generalization 

2. Ordinary complex numbers. Let a and b be any real numbers. 
Then the complex number a + bi is said to have the norm a 2 + b 2 . For- 
mula (1) evidently expresses the property that the norm of the product 
r + si of the complex numbers a + bi and a + pi equals the product of 
their norms. 

To prepare the way for our introduction to quaternions and Cayley's 
algebra, we shall present briefly W. R. Hamilton's definition of complex 
numbers by means of couples of real numbers. Two couples (a, 6) and 
(a, 0) are called equal if and only if a = a, b = /?. Addition and multi- 
plication are defined by 

(a, b) + (a, 0) = (a + a, b + 0), (a, &)(«, 0) = (r, s), 

where r and s are given by (1). If to is any real number, we define 
m(a, b) and (a, b)m to be (ma, mb). Writing 1 for (1, 0) and i for (0, 1), 
we have 

(a, b) = (a, 0) + (0, b) = o(l, 0) + 6(0, 1) = a + bi. 

The previous definition of addition and multiplication of couples gives 

(a + bi) + (a + /Si) = a + a + (b + p)i, (a + bi)(a + pi) = r + si. 

155 
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3. Quaternions. Consider quadruples (a, b, c, d) of real or complex 
numbers a, b, c, d. Define addition and multiplication by 

(a, b, c, d) + (a, P,y,S) = (a + a,b + P, c + y, d + 5), 

(a, b, c, d) X (a, p, y, 5) = (A, B, C, D), 
where 

A = aa — bp — cy — dS, B = ap -f- ba + c8 — dy, 



(2) 



C = ay - bS + ca + dp, D = aS + by - cP + da. 



No attempt will be made here to explain why we select these values for 
A, • • •, D; it is not our purpose to explain how quaternions were discov- 
ered or how they may be made to enter naturally,* as we aim merely to 
give a logical basis for quaternions. Consider the four particular quad- 
ruples 

1 = (1, 0, 0, 0), i = (0, 1, 0, 0), j = (0, 0, 1, 0), k = (0, 0, 0, 1), 

called the units. Define m(a, b, c, d) or (a, b, c, d)m to be (ma, mb, mc, 
md), where m is any complex number. Then 

(a, b, c, d) = (o, 0, 0, 0) + • • • + (0, 0,0, d) = a + bi + cj + dk, 

i 2 = j 2 = k 2 = — 1, ij = k, ji = — k, 

(3) 

jk = i, kj = - i, ki = j, ik = - j. 

Henceforth we discard the quadruple notation and employ 

q = a + bi + cj + dk, Q = a + pi + yj + dk, 

called quaternions. In view of our earlier definitions, their sum is a + a 
+ (b + p)i + • • • and their product is A + Bi + Cj + Dk, where A, 
■ • ■ , D have the values (2) . This product may be found by performing 
the multiplication as in formal algebra, care being taken not to permute 
two factors i,j,k, and then simplifying the result by use of (3). For 
example, (i + 2j)(j + k) = k — j — 2 + 2i. Note that, while multipli- 
cation is not commutative, it is associative since (ij)k = — 1 = i(jk), etc. 
The quaternion q' = a — bi — cj — dk is called the conjugate to q. 
We readily verify that qq' = q'q = a 2 + b 2 + c 2 + d 2 , which is called the 
norm N(q) of q. For the moment, let a, b, c, d be real numbers, so that 
q is a real quaternion; if q 4= 0, then N(q) # 0, and q has the inverse 
g -1 = q'/N(q). Thus, if g # 0, qQ = qi has the unique solution Q = q~^qu 
and Qq = q^ has the unique solution Q = qiq* 1 , so that both right-hand 

* This topic is presented in an elementary manner in Dickson's Linear Algebras, Cambridge 
University Tract No. 16, pp. 9-12, and from another standpoint in his article "On the relation 
between linear algebras and continuous groups," Bull. Amer. Math. Soc, 22, 1915, 53-61. 
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and left-hand division are always uniquely possible if the divisor is a real 
quaternion not zero. 

The conjugate of qQ equals the product Q'q' of the conjugates of the 
factors taken in reverse order, as shown by interchanging the Roman and 
Greek letters in the sums (2) and afterwards changing the signs of b, c, 
d, P, y, 5. 

The norm of qQ is qQ • Q'q' by definition. By the associative law, this 
may be written q{QQ')q'. Since QQ' is an ordinary number, it is commu- 
tative with q' in view of our earlier definition of m(a, b, c, d) and (a, b, c, d)m. 
The result is now the product of the norms qq' and QQ' of q and Q. Hence 
the norm of a product of two quaternions equals the product of their 
norms, i. e., 

(a 2 + b 2 + c 2 + d 2 )(« 2 + /3 2 + t 2 + 8 2 ) 

(4) 

= A 2 + B 2 + C 2 + D 2 (A, ■ ■ •, D as in (2)). 

Much earlier than Hamilton's invention of quaternions in 1843, Euler* 
discovered formula (4) while investigating the elegant theorem that every 
positive integer is a sum of four integral squares, the theorem following 
from (4) if proved for every prime number; he also used (4) in his later 
paper on orthogonal substitutions. 

4. Cayley's algebra. A. Cayleyf defined an algebra with the 8 units 
1, i u • ■ -, i 7 , such that i\ = — 1, • • •, i 7 2 = — 1, 

HH = U = — tji'i, iftz = i\ = — i&i, Hi\ = ii = — iiiz, 

and six similar sets of six relations with 1, 2, 3 replaced by 1, 4, 5; 6, 2, 4; 
6, 5, 3; 7, 2, 5; 7, 3, 4; 1, 7, 6; respectively. Then 

(z + Ziii H \- x 7 i 7 )(x ' + xi'ii -\ f- x 7 'i 7 ) = A + Aii t + (- A 7 h, 

where, if we employ the abbreviations jk = XjX k ' — x k x/, Oj = x x/ + XjX ', 

A = x Xo' - XiXi - ... - aW, At = 23 + 45 + 76 + 01, 

A 2 = 31 + 46 + 57 + 02, A s = 12 + 65 + 47 + 03, 

A 4 = 51 + 62 + 73 + 04, A 5 = 14 + 36 + 72 + 05, 

A 6 = 24 + 53 + 17 + 06, A 7 = 25 + 34 + 61 + 07. 

He called ~2Xi 2 the modulus (norm) of x + ■ ■ ■ + x 7 i 7 and stated that the 

* Corresp. Math. Phys. (ed., P. H. Fuss), I, 1843, 452, letter to Goldbach, May 4, 1748. 
Novi Coram. Acad. Petrop., 5, 1754-5, 3; 15, 1770, 75; Coram. Arith! Coll., I, 230, 427. 

fPhil. Mag., London, (3), 26, 1845, 210 [30, 1847, 257-8]; Coll. Math. Papers, I, 127 [301]. 
In Ai his 87 is a misprint for 73. 
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norm of a product equals the product of the norms of the two factors: 

(5) (£xA(£x i ' 2 ) = j:A i *. 

\ i=u / \ i=0 / i=0 

The last result, as well as another important property of the algebra, 
can be proved without computation by representing the algebra as a 
quasi-binary algebra.* Since 1, i u i 2 , is satisfy the relations (3) for the 
quaternion units, we may replace them by 1, i, j, k. Then the remaining 
four units are e = u, ie = i & , je = i 6 , ke = i 7 . Hence every number of 
the algebra is of the form q + Qe, where q and Q are linear functions of 
1, i,j, k and hence are quaternions. It can be verified that Cayley's 49 
relations, giving the product of two equal or distinct units i x , • • • , i 7 , are 
together equivalent to the single formula 

(6) (q + Qe)(r + Re) = qr - R'Q + (Rq + Qr')e, 

where r' and R' are the quaternions conjugate to r and R. The reader 
need not verify the equivalence stated, but may take (6) as the rule of 
multiplication for the numbers of the algebra to be considered henceforth, 
since Cayley's algebra has been introduced here merely for historical back- 
ground and will not be further employed in his form. 

Define the norm of q + Qe to be qq' + QQ', which is a sum of 8 squares. 
Taking r = q', R = — Q, in (6), we get 

(q + Qe)(q' - Qe) = qq' + QQ', 

so that the norm of q + Qe is its product by its conjugate q' — Qe. Since 
multiplication does not here obey the associative law, we cannot conclude 
at once, as we did for quaternions in §3, that the norm of a product equals 
the product of the norms of the two factors. However, we obtain a short 
proof by use of a device. Express the right member of (6) in the form 
t + Te by setting 

(7) t = qr - R'Q, T = Rq + Qr'. 

Its norm tt' + TT' is seen, by direct multiplication and use of the fact 
that the norm of qr is r'q', to equal a — |3 + y, where 

a = RqrQ' + Qr'q'R', = qrQ'R + R'Qr'q', 

y = qrr'q' + R'QQ'R + Rqq'R' + Qr'rQ' = (qq' + QQ')(rr' + RR'). 

The last equality is a consequence of the fact that rr' is an ordinary num- 
ber and hence can be interchanged with q', etc. Our device occurs in 
the proof that a = j3. Note that the conjugate of the first term of a 

* Dickson, Trans. Amer. Math. Soc, 13, 1912, 72; Linear Algebras, 1914, 15. 
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equals the second term of a, so that a is an ordinary number and hence is 
commutative with every quaternion. Hence a = R'aR -j- RR', which is 
seen to equal /3. In the excluded case R = 0, evidently a = /3 = 0. 
Hence the norm of the product (6) equals the product of the norms of the 
factors. Thus we can write down an 8 square formula of type (5). 

Moreover, both right-hand and left-hand division except by zero is 
always possible and unique in our algebra composed of the numbers q 
+ Qe, provided we restrict q and Q to be real quaternions. Of the two 
types of division consider that in which the second factor r -f- Re and the 
product t + Te are given, while the first factor q + Qe is to be found. 
Thus we seek to solve equations (7) for q and Q. Multiply the second 
equation (7) by r on the right and replace qr by its value from the first 
equation; we get 

(rr' + RR')Q = Tr - Rt. 

Again, multiply the first equation byr' on the right and eliminate Qr'; thus 

(rr' + RR')q = tr' + R'T. 

Since rr' + RR' equals the sum of the squares of eight real numbers, it 
is zero if and only if r = R = 0. Similarly, equations (7) can be solved 
for r + Re unless q = Q = 0. 

We have now accomplished one of the aims of the paper, having ex- 
hibited linear algebras in 2, 4 and 8 units for which the norm of a product 
equals the product of the norms of the factors (thus giving the 2, 4 and 8 
square theorems), and such that, if the coordinates of the numbers of 
the algebra be restricted to be real numbers, both right-hand and left- 
hand division except by zero are possible and unique. While the three 
algebras have in common these two fundamental properties, they differ 
in other respects. For complex numbers multiplication is both commu- 
tative and associative, for quaternions it is associative but not commuta- 
tive, for Cayley's algebra of 8 units it is neither commutative nor associa- 
tive. What additional properties must be given up to obtain a similar 
linear algebra in more than 8 units? We shall prove in §5 that there 
exists no linear algebra in more than 8 units for which the norm is a sum 
of squares and the norm of a product equals the product of the norms of 
the factors. 

5. Hurwitz's Theorem.* We seek the values of n for which there exists 



* Gottingen Nachrichten, 1898, 309-316. Since experience shows that graduate students 
fail to follow various steps merely outlined by Hurwitz, we shall here give the proof in detailed, 
amplified form. As we shall employ (a,-,-) to denote a matrix and not a linear transformation, we 
must invert the order of factors in his products. 
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an identity (as to the x's and y's) of the form 

(8) (Xj 2 + • • • + X„ 2 )(?/l 2 + • • • + Vn 2 ) = Zl 2 + • • • + Zn\ 

where z\, • • ■, z„ are linear in x u • ■ •, x n and also in y u • ■ •, y n . Let 

(9) Zi = a,i2/i + • • • + a in y n (i = 1, ■ ■ •, n), 
where the ay are linear functions of xi, ■ ■ ■ , x n . We employ the matrices 

(On O12 • • • o x „ \ /On o 2 i • • 

a 2 l O22 • ■ ■ 2n J A' = I ffll2 ° 22 ' ' 

Orel O n 2 - - - Orere / \"l« 02n " * 

where A' is derived from A by the interchange of its rows and columns, 
and is called the conjugate (or transposed) of A. In case the diagonal 
elements ay all equal a and the elements not in the diagonal are all zero, 
we shall write al for A, where / is the unit (or identity) matrix and has 
the property that IB = BI = B for every matrix B of n rows and n 
columns. The quadratic form 

n 

^2 bijZiZj = buZx 2 + 2&12Z1Z2 + &22Z2 2 + • • • (by = bji) 

is said to have the matrix B = (by), whose ith row is b a , b i2 , • ■ • , b in . 
If we replace the variables z h • • ■ ,z n by the expressions (9), we evidently 
obtain a new quadratic form in the variables y lt •••,y n ', its matrix is 
known* to equal A'BA. In particular, let the quadratic form be Z1 2 
+ • • • + z„ 2 , whose matrix is B = I; then the quadratic form derived 
by replacing z x , • ■ •, z„ by the expressions (9) has the matrix A' A, a fact 
which can be verified at once without making use of the standard theorem 
just quoted. Now we desire that the resulting quadratic form in y u • • ■ , 
y n shall be the left member of (8), whose matrix is al, where a = X1 2 
+ • • • + x„ 2 . Hence there exists an identity (8) if and only if there exist 
n 2 linear functions ay of Xi, • • •, x„, whose matrix (ay) is denoted by A, 
such that 

(10) A' A = (Xx 2 + • • • + x„ 2 )7. 

Since each element of matrix A is a linear function of x u • • • , x„, and 
since the sum of several matrices is a matrix whose elements are the sums 
of the corresponding elements in the matrices added, it follows that A 
= X1A1 + • • • + x„A n , where A u ■ ■ •, A n are matrices with constant ele- 
ments. Thus in A' A the coefficient of x„ 2 is A n 'A n , which equals / by 
(10). Let B { = A n 'Ai (i = 1, • • •, n — 1), whence A t = A„By A/ 

* Bficher, Introduction to Higher Algebra, p. 129. 
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= B/An, and A 'A equals 

(XiBi + ■ ■ • + X„_iB;,_ 1 + X n )A n ' • A n (XiBi + • • • + X n -iB n -i + X n ). 

Since A„'A n = I, (10) becomes 

(11) (xiBi' + ■■■ + x„_iS;_, + x^ixxBx + ■ • ■ + x n -iB n -i + x n ) 

= (Xi 2 + • • • + X n *)I. 

Thus Bi'Bi = I, Bi' + Bi = 0, B/B k + B k 'Bi = 0, whence 

(12) Bi' Bi, Bt' = - I, BiB k = - B k Bi 

(i, k = 1, • ■ ■, n — 1; i =j= k). 

A matrix B = (& y ) is called symmetric if b^ = by, and skew symmetric 
if bn = — by for every i, j; thus B is symmetric if and only if B' = B, 
and skew-symmetric if and only if B' = — B. The latter condition im- 
plies that b = ( — ) n b if b is the determinant of the matrix B of n rows and 
n columns. Thus b = if n is odd. By the first two equations (12), 
Bi is skew-symmetric and its determinant is not zero, so that n is not odd. 
Hence there exists no identity (8) if n is odd. In what follows, we assume 
that n is even. 

Our next step is to prove that at least half of the matrices 

(13) /, B H , BiJ5i V BifiiJB^, • • •, B1B2 ■ ■ ■ B n -i 

(t'i < n, i'i < it < n, • • •) 

are linearly independent. There are 2" _1 such products since any one 
product either contains Bi or does not, • • •, and either contains -B n -i or 
does not. Let G = B it • • • B ir be one of the matrices (13) ; it is symmetric 
if r = or 3 (mod 4), and skew-symmetric if r = 1 or 2 (mod 4), since 
by (12) 

G' = B^ • • • Bit = (- \YB ir ...B h = (- 1)'G, 

where s = r + r — 1+r — 2 + • • • + 1 = r(r + l)/2 is even if r = 0, 3 
(mod 4), but odd if r = 1, 2 (mod 4). In particular, a product of two dis- 
tinct B's is skew-symmetric. 

Consider the possible linear relations (with constant coefficients not 
all zero) which hold between the matrices (13). Such a relation R = 
is called irreducible if it is not possible to express R in the form R = Ri 
+ Rt, where Ri = and R 2 = represent two linear relations holding 
between our matrices such that no one of these matrices (13) occurs as a 
term of both Ri and R 2 . In particular, an irreducible linear relation does 
not involve both symmetric and skew-symmetric matrices, since it could 
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then be written in the form M = S, where M is the aggregate of the sym- 
metric matrices and £ is the aggregate of the skew symmetric matrices^ 
whence M' = S', M' = M, S' = - S, giving M = 0, S = 0. 

Let R = be any irreducible linear relation between the matrices (13). 
By multiplying R by the product of a constant and a suitably chosen 
matrix (13), we get a new linear relation p = 0, one term of which is I 
and all the remaining terms are products of matrices (13) by constants. 
Thus if ABzBz is one term of R, we use the multiplier — \B 2 B 3 . We 
need also to know that if we multiply the matrices (13) on the left by any 
one (say M) of them, the products form a permutation of those matrices 
each prefixed with the factor + 1 or — 1. This is evident when the multi- 
plier is Bi, since the product will contain or lack Bi according as the multi- 
plicand (13) lacks or contains B u in view of Bi 2 = — I. If the multi- 
plier is B 2 , we first replace BiB 2 • • • by — B 2 Bi • ■ • and see that the for- 
mer argument applies. After proving in this manner our statement when 
the multiplier is any B { , we see that it holds when the multiplier is any 
product of the B's. Returning to our new relation p = 0, we note that 
it also is irreducible, since by multiplying it by a product of a constant 
and a suitable matrix (13) we recover our initial relation R = 0, which 
was assumed irreducible. Hence p = is an irreducible relation 

J- = SCjij^-Di^-jij + Z/di 1 i i i z i l Bi 1 Bi i BiJSi i + • • • 

involving exclusively symmetric matrices (13), so that no term contains a 
single Bi or a product of only two B's. Multiply all the terms of our re- 
lation by Bi on the right; we obtain an irreducible relation which there- 
fore involves only skew-symmetric matrices (13), one term being Bi. 
Since a product of four distinct B's is symmetric, we conclude that c,-,^, 
is zero if i is distinct from i u i 2 , i$. Since i may have any value Si n — 1, 
we have c = unless 3 = n — 1. To prove that every d = 0, take 
i = i 4 ; then the coefficient of — dB^B^B^ is zero. The method used to 
prove c = applies when the number r of factors B is = 3 (mod 4) and 
r < n — 1, since r + 1 = 0. The method used to prove d = applies 
when r = (mod 4), since r — 1 = 3. Hence if our relation exists, it 
has the form 

J = kBiB 2 • ■ ■ B_i. 

Since each member is a symmetric matrix, n — 1 = or 3 (mod 4) . But 
n is even. Hence n = (mod 4). As in the discussion of G, below (13), 
the square of Bi • • ■ B r is (— 1)"7, where s = r(r + l)/2. Hence k 2 = 1. 
Thus the 2 n ~ 1 matrices (13) are linearly independent if n = 2 (mod 4); 
while for n = (mod 4) they are either linearly independent or are connected 
by the relations which arise from I = ± BiB 2 • • ■ -B«-i by multiplication by 
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the various matrices (13), but are connected by no further irreducible linear 
relations. 

To illustrate this result, let n = 4. Then the 8 matrices 

I, B\, Bi, Bi, B1B2, BiBs, B2B3, BiB 2 B s 

are either linearly independent or are connected by only four irreducible 
linear relations; 

J = ± B1B2B3, Bi = =F B2B3, B2 = ± BiBs, Bi = =F B1B2. 

The latter express BiB 2 Bz, B 2 B 3 , B1B3, BiB?. linearly in terms of I, B lt 
B 2 , B it which are therefore in all cases linearly independent. 

For any n, one of the reduced products of I and B t ■ • ■ B„_i by any 
matrix (13) evidently contains fewer than half of the B's and the other 
contains more than half of the B's. Hence if irreducible linear relations 
exist, they serve merely to express the latter products in terms of the for- 
mer. Thus in every case, the 2" -2 matrices (13) which are products of 
at most (n — 2)/2 factors B are linearly independent. 

But if we are given any n 2 + 1 matrices (aij ik) ) each with n rows and 
n columns, we can find numbers x* not all zero such that 

E **(a«,<») = 0, 
i. e., 

7*2 + 1 

E x*a»/*° = (i, j = 1, • • •, n), 

k=\ 

since n 2 linear homogeneous equations in n 2 + 1 unknowns x k have solu- 
tions not all zero. 

Hence 2"~ 2 ^ n 2 . This is satisfied if n S 8, but fails if n = 10. But 
if it fails for n = m, it fails for n = m + 1, since 

2 m+i-2 = 2.2 m - 2 > 2to 2 > (w + l) 2 

if (to — l) 2 > 2, and hence if to S 3. We have now proved that n ^ 8. 
The case n = 6 is readily excluded. Then the 2 5 matrices (13) are 
linearly independent. But 5 + 10 + 1 of them are skew-symmetric 
(those with 1, 2 or 5 factors B). Between any 16 skew-symmetric six- 
rowed square matrices there exists a linear relation : 

16 16 

E **(&«<») = 0; X)x t &«<» = (i,j = 1, • • •, 6; i < j), 

k=l k=\ 

it being now necessary to examine only the 15 terms to the right of the 
main diagonal. But 15 linear homogeneous equations in 16 unknowns 
Xk have solutions not all zero. 
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Theorem. Except for n = 1, 2, 4, 8, there exists no identity (8) ex- 
pressing the product (xi 2 + • • • + x n 2 )(yi 2 + • • • + y n 2 ) as a sum of the 
squares of n bilinear functions of X\, • • •, x n and y i} ■ ■ ■, y n . 

History of the Subject. 

6. Gauss* remarked that the four square formula (4) is expressed in 
a simple way by 

(Nl + Nm)(N\ + Nix) = N(l\ + mix) + NQix' - mk'), 

where I, m, X, n, X', \x are complex numbers, X' being conjugate to X, and 
ix to ix, while Nl denotes the norm of I (§2). 

7. C. F. Degenf extended Euler's formula (4) to eight squares : 

(P 2 + Q 2 + W + S 2 + T 2 + U 2 + V 2 + X 2 ) 

X (p 2 + q 2 + r 2 + s 2 + t 2 + u 2 + v 2 + x 2 ) 

= (Pp + Qq + Rr + Ss + Tt + Uu + Vv + Xx) 2 

+ (Pq - Qp + Rs - Sr + Tu - Ut + Vx - Xv) 2 

+ (Pr - Qs - Rp + Sq=F Tv ± Ux ± Vt ^ Xu) 2 

+ (Ps + Qr - Rq - Sp ±Tx±Uv^Vu^ Xt) 2 

+ (Pt - Qu ± Rv T Sx - Tp + Uq T Vr ± Xs) 2 

+ (Pu + Qt T Rx =F -Sr - Tq - Up ± Vs ± Xr) 2 

+ (Pv - Qx =F Rt ± Su ± Tr T Us - Vp + Xq) 2 

+ (Px + Qv±Ru±St^ Ts=F Ur - Vq - Xp) 2 . 

He stated [erroneously as we saw in §5] that there is a like formula for 
2" squares. For the case of 16 squares he gave the literal parts of the 
16 bilinear functions, but left most of the signs undetermined, saying that 
the only difficulty is the prolixity of the ambiguities of signs. This paper 
has been overlooked by all subsequent writers on the subject. 

8. J. T. Graves} communicated to W. R. Hamilton Jan. 18, 1844 
(correcting some errors in signs in the formula communicated Dec. 26, 
1843), a formula which differs from Cayley's (5) only in the interchange of 
6 and 7, and a second formula which becomes Cayley's on writing x , 
• • • , x 7 for a, b, ■ • • , h. Hence Graves's formulas need not be inserted 

* Posthumous MS., Werke, 3, 1876, 383-4. 

t Mem. Acad. Sc. St. Petersbourg, 8, annees 1817-8 (1822), 207-219. There is a misprint in 
the sign of his term ± Rt, here corrected. 

JProc. Roy. Irish Acad., 3, 1845-7, 527-9; Trans. Roy. Irish Acad., 21, II, 1848, 338-341; 
Phil. Mag. London, (3), 26, 1845, 320. 
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here. At first he expected that it would be possible to give an extension 
to 2" squares. 

9. J. R. Young's* formula, with s, t, u, v, y, z, w, x replaced by a h • • • , 
a 8 , is 

(14) (ga^)^ 1 ) = (-Earn) 2 

+ (12 + 34 + 56 + 78) 2 + (13 + 42 + 57 + 86) 2 

+ (41 + 32 + 58 + 67) 2 + (15 + 62 + 73 + 48) 2 

+ (16 + 25 + 38 + 47) 2 + (17 + 82 + 35 + 64) 2 

+ (18 + 27 + 63 + 54) 2 , 

where ij denotes atuj — a.i<ij. It was admitted to be equivalent to Graves's 
formulas. Youngf stated that a like formula holds for 2" squares, but 
soon afterwards admitted that this is erroneous, saying that he was pre- 
pared to prove that the proposition does not hold beyond 8 squares. 

Young} gave a long discussion to show that the extension to a sum 
$i6 of 16 squares is impossible. He exhibited a special relation &-&&W 
= Su" in which the roots of 8 of the squares in S u are proportional to the 
roots of 8 of the squares in Su'. He§ noted that, for k = 2, 4, 8, a product 
of a sum of km squares by a sum of kn squares can be expressed as a sum 
of kmn squares. 

10. Cayleyjl investigated the possibility of a formula for 2" squares by 
introducing 2" — 1 symbols a , &o, ■ ■ • , not assumed to be commutative, 
but such that a 2 = b 2 = • • • = — 1 and 

boCo = ± a = — Co&o, c ao = ± b = — aoC , a bo = ± Co = — &o«o. 

Denoting this set of six equations by a &oCo = ±, let also a doCo = ±, 
etc., where the sign is not necessarily the same as before, while the sys- 
tem of triples contains each duad once and but once, and the signs are to 
be chosen at will. Then 

(w + aa + bb + • ■ -)(wi + a,ia + 6i6 +•••) = w 2 + a 2 a + b 2 b + ■ • •, 

where w 2 , a,z, • • ■ are linear and homogeneous in w, o, • • • and in w u a u 

• • •. Assume (I) that if any two triples with a common element, e aob 

and e Codo, occur in the system, there occur also foa c 0) fod o b , gododo, goboCo ; 

* Proc. Roy. Irish Acad., 3, 1845-7, 526-7. 
t Phil. Mag., London, (3), 30, 1847, 424-5; 31, 1847, 123. 

t Trans. Roy. Irish Acad., 21, II, 1848, 311-338. Outline in Proc. Roy. Irish Acad., 4, 1847 
-50, 19-20. 

§ Phil. Mag., London, (3), 34, 1849, 114. 

|| Phil. Mag., London, (4), 4, 1852, 515-9; Coll. Math. Papers, II, 49-52. 
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(II) that for any two pairs of triples, such as e ao&o, e c do and/ aoCo, fodob 0) 
the products of the signs of the triples in the first pair is the same as that 
in the second pair. Then 

O 2 + a 2 + b 2 + ■ ■ -)(wi 2 + «i 2 + V H ) = «>2 2 + a 2 2 + & 2 2 + • • •■ 

The converse was not proved, but it was stated that conditions (I) and 
(II) afford a complete test for the possibility of the 2 n square theorem. 

T. P. Kirkman,* to whom Cayley had communicated privately the 
preceding test, verified that, for 15 elements a,b, • • • , triples can be chosen 
so that (I) is satisfied, but that (II) then involves a contradiction. 

11. F. Brioschif showed that, if re is even, the square of the determi- 
nant A = | a,,- 1 of order re is a skew-symmetric determinant L = | Zy | of 
order re with the general element 

In = tt r iffl«2 — tt r 2<Isl T dr3<J>si — 0>riO>a3 T" " " ' ~\~ dm-llsn — dmO-sn—l = — hr- 

Similarly, the square of C = | c,-,- 1 is | py |, where p rs = c r ic s2 — • • •. Let 

n 
AG = [ Ay |, A rs = / , ClriCsi, I Ay I =| Lij I, L rs = A r lA B 2 — • • •. 

If the a's and c's are such that 

(15) hi = I34 = • • • = ln-ln = t, P12 = • • • = Pn-lr, = W, 

while the remaining Uj and py are zero, it is proved that L i2 =i L 34 = • • • 
= L n _i„ = tu, and that the remaining La are zero. Now let re = 8 and 
take a ti = an, ay = — ay (i + j) except for a 1& = a 5 i = a 26 = a 62 = 037 
= a 73 = a 48 = a 84 and take also 

ai 2 = &43 = #56 = tt87, <^13 = C 2 4 = Os7 = 068, ffl 14 = a 32 = Oss = 076, 

ai6 = a 4 7 = «5 2 = a83, an = a 2 s = a§3 = a64, ai$ = a 3 6 = as 4 = 072. 

Assume like relations between the d,. It is stated erroneously that rela- 
tions (15) and the analogous relations between the Ay hold, so that 

SAy 2 = tU, t = Say 2 , M = SCy 2 (j = 1, • • •, 8). 

Although Z12 = hi = Z 56 = hs = Say 2 , it was pointed out by E. Sadunf 
that 

Jie = 2(anai5 + ai 2 aie + a 13 ai 7 + ai 4 ai 8 ) 4= 0, 

so that we cannot make t = Say 2 . In a footnote, Sadun reconstructed 
Brioschi's proof, and obtained (14) with 5 and 7, 6 and 8 interchanged. 

* Phil. Mag., London, (3), 33, 1848, 447^159, 494-509; (3), 37, 1850, 292-301. 
t Jour, fur Math., 52, 1856, 133-141; Opere Mat., V, p. 511. 
t Periodico di Mat., 14, 1899, 125-139; and pamphlet of 1896. 
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12. A. Lebesgue* gave an 8 square formula, communicated to him 
to Prouhet, which apart from signs becomes Cayley's formula (5) if we 
write x , • ■ ■, x 7 for a,b, • • • , h. 

13. A. Genocchif concluded that sums of 2" squares repeat under 
multiplication by an erroneous argument (false even for n = 2) based 
upon sums of two squares. The error was pointed out by Sadun (§11) 
and earlier by A. Puchta,J who interpreted the correct 8 square formula 
by means of regular bodies with 9 vertices in space of 8 dimensions. 

14. E. Mathieu§ expressed Euler's identity (4) in the form 

SSi = S/', Si = x 2 + x x 2 + x a 2 + x 2 l+tc , w 2 + w + 1 = (mod 2), 

while x w " and x'{ +w are derived from x" by the substitution (z, wz), viz., 
(1, w, 1 + w), on the subscripts. But Si" is unaltered also by (0, w, 
1 + w). Hence of the 24 permutations on the four subscripts, 12 give 
one decomposition into 4 squares, and 12 give another. 

For w 3 + w + 1 = (mod 2), Cayley's formula (5) can be expressed 
in the form 

SsS*' = S t ", S S = SZy 2 , X " = HXjX/, 

\ Xl^ W 2X-u;2 X W +w2Xl+ w + w 2 "T~ Xl-\- W ^- W 2X w ^- W i, 

where j ranges over the eight values 0, • • • , 1 + w 2 appearing in 

s = (0)(1, w, w 2 , 1 + w, w + w 2 , 1 + w + w 2 , 1 + w 2 ). 

The remaining x" are derived from X\" by applying this substitution 
s, which may be written in the form (w z , w z+1 ), the signs of the terms of 
x/' being determined so that the terms occurring in the above *S 4 " occur 
with the same signs in S$". Now Xi" is unaltered by (w z , w 2z ), while 
Zo" 2 , • • •, x 7 " 2 are permuted by (w z , w llz ), where 1/z is replaced by the in- 
teger congruent to it modulo 7. Hence any symmetric function of these 
8 squares is unaltered by the 3-7-8 substitutions 

(w z , w*'), z' = ^^r d , ad-bc = 1, 2, 4 (mod 7). 

It is stated that these results cannot be extended to more than 8 
squares. 

* Exercices d'analyse num6rique, 1859, 104; Introduction a la thcorie des nombres, 1862, 65. 
t Annali di Mat., 3, 1860, 202-5; Giornale di Mat., 2, 1864, 47-48. 
t Sitzungsber. Ak. Wiss. Wien (Math.), 96, II, 1887, 110-133. 
§ Jour, fur Math., 60, 1862, 351-6. 
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15. J. J. Thomson* verified Young's formula (14) by means of rela- 
tions like 

(16) 12-34+ 13-42 + 41-32 = 0. 

16. E. Lucas statedf that there is a relation between the formula 
expressing the product of two sums of n squares as a sum of n squares for 
n = 4, 8, 16, etc., and Sylvester's| square diagram formed of an equal num- 
ber of white cases and black cases, such that for any two lines or two col- 
umns the number of variations of colors is always equal to the number of 
permanences. If, in the accompanying diagram, we replace each white 
case by a plus sign and each black case by a minus sign, we are led to 
Euler's formula (4). 



17. S. Roberts§ argued that a 16 square formula is impossible. He 
assumed in effect that an m square formula must be of the type 

(17) 2> 2 ( 2> 2 ) = Z (oiCfl + • • • + a m c im Y, 

where Cn, • • -, c im and c u , • • -, c mi are permutations of ± c u • • -, ± c m , 
and that the formula reduces to a \m square formula by setting a { = c { 
= (i > to/2). In building the to square formula from the \m square 
formula, he made free choice between the letters not already in the scheme. 
He derived an unique formula for m — 4 and for m = 8, but found after 
a tedious examination a contradiction for to = 16. 

18. Cayleyj| considered the linear algebra with the units E = 1, E u 
■ • •, Et, where E { 2 = - 1 (i = 1, • • -, 7) and 

E1E2E3 = «!, E\E^Ef, = « 2 , EzEiEb = €4, EzE^Ei = « 6 , 

E\E(,Ei = €3, EiE^Ei = 65, EgEsEs = ey, 

* Messenger Math., 7, 1877-8, 73-74. 

f Assoc, franc, av. so., 6, 1877, 213-4. 

t Math. Quest. Educ. Times, 10, 1868, 74-6, 112 (diagrams for 8 X 8 squares and 16 X 16 
squares). Cf. M. Jenkins, ibid., 14, 1871, 22-25. 

§ Quar. Jour. Math., 16, 1879, 159-170. 

|| Amer. Jour. Math., 4, 1881, 293-6; Coll. Math. Papers, XI, 368-371. Incomplete summary 
in Johns Hopkins University Circulars, 1882, 203. 
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each €,- being 1 or — 1, and the first symbol denotes the six equations 

E\Ei = «ix?3 = — E2E1, E<iE% = t\E\ = — E 3 Ei, E$E\ = e%E 2 = — E\E$. 

For no values of the e's is the algebra associative. We may set 

(SaiEiXSa/Ei) = Xa/'Et (i = 0,1, • ■ ■, 7). 

Without loss of generality we may take ei = e 2 = «3 = + 1. Then 

(Za,- 2 )(Za/ 2 ) = Sa/' 2 , 

if and only if — e 4 = e 5 = « 6 = «7- In such an algebra, E x E 2 • E 3 
= Ei-E 2 E 3 and similarly for each of the seven triads above. For the re- 
maining 28 triads, EiE r E k = — E { -EjE k . [We pass from one of these 
two algebras to the other by changing the signs of E 2) • • • , E 7 . If we take 
€4 = — 1 and change the sign of E 7 , we get Cayley's earlier algebra (§4).] 

19. Cayley* remarked that (16) establishes Euler's identity! 

(Za^XZyi 2 ) - (2z,-2/,) 2 = (12 + 34) 2 + (13 - 24) 2 + (14 + 23) 2 . 

The first step in forming this identity is to arrange the duads into a syn- 
thematic form: 12-34, 13-24, 14-23. The next step is to determine the 
signs. For 8 elements there is a single such synthematic arrangement; 
if 34, 56, 78 and each lj are taken with positive signs, only one sign remains 
arbitrary, so that there are only two final schemes. For 16 elements, we 
have first to form 15 lines each containing the numbers 1, • • -,16 in 8 
duads, no duad being repeated. Only four types are found; for each it 
is found to be impossible to choose the signs. Cayley states that earlier 
writers had tacitly assumed that only one of the four types is possible and 
hence had not given a complete proof of the non-existence of the 16 square 
theorem. The question of the distinctness of the four types, apart from 
notation, was mentioned, but not discussed, by Cayley. 

20. S. Robertsf remarked that Cayley's four types are all equivalent. 
But his directions for deriving the first from the second type are incorrect. 
Besides interchanging 13 with 14, and 15 with 16, and interchanging col- 
umns 13 and 14, and rows 15 and 16, it is necessary to interchange also 
rows 13 and 14, and columns 15 and 16. He indicated how his own pro- 
cess can be used to produce the four (equivalent) types. 

21. F. Studnicka§ employed the product of two determinants: 

a b 
-V a' 





x y 




ax + by 


— ay' + bx' 




-y' x' 




a'y — b'x 


a'x' + b'y' 



* Quar. Jour. Math., 17, 1881, 258-276; Coll. Math. Papers, XI, 294-313. 

f Quoted in Math. Quest. Educ. Times, 75, 1901, 40. 

t Quar. Jour. Math., 17, 1881, 276-280. 

§ Sitzungsberichte K. Bohm Gesell. Wiss. Prag, 1883, 475-481. 
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Taking a' = a, • • • ,y' = y, we get (1). Next, let a' be the conjugate to 
the complex number a, ••■,y' the conjugate to y; we get Euler's (4). 
But he erred in employing the same formula when a' and a are conjugate 
quaternions, • • •, y' and y conjugate quaternions, to deduce the 8 square 
theorem, since he overlooked the fact that the initial formula holds only 
when multiplication is commutative [Vahlen, §27]. 

22. X. Antomari* wrote (a,ft-) for aft, — piaj and employed the iden- 
tity 

D = (ZdiXi) (2&4/ f ) - (Sa,-2/i)(S6ia;,) = ~Z{a t b J ){x i y J ) (i,j = 1, • • -, 4; j > i). 

In view of (16), written in a, b and again in x, y, we get 

D = { (a!& 2 ) + (x 3 2/ 4 ) ) { O12/2) + (a 3 & 4 ) } 

+ { O2&3) + O12/4) } { O22/3) + (0164) } + { O1&3) + O42/2) } { O12/3) + (a 4 6 2 ) } • 

Taking aj and x, to be conjugate complex numbers and also bj and y,, for 
j = 1, • • •, 4, we get an 8 square formula. 

23. E. Lucasf stated that the determinant of the 8 equations 

ax + by + cz + dt + ep + fq + gr + hs = X, ■ ■ ■ , — hx + ■ ■ • + as = S 

is the fourth power of A = a 2 + • • • + h 2 . To solve the equations, mul- 
tiply them by a, • • • , h, taken with proper signs, and add. We get 

Ax = aX - bY - cZ - dT - eP - fQ - gR - hS, ■ ■ ■, 

As = hX + ■ ■ ■ + aS. 
Squaring these and adding, we get Sa 2 -Sx 2 = SX 2 . 

24. G. ArnouxJ argued the impossibility of a 2" square formula for 
n > 3. 

25. Teilhet, de Montessus and Boutin § gave special numerical ex- 
amples of S r S r ' = S T " for r = 16 and r = 32, where »S r is a sum of r 
squares. 

26. E. Sadun|| discussed m square formulas of the type (17). Since 
the product terms in a p a t must cancel if p + t, we have 

Without loss of generality we may assume that the first row and first col- 
umn of the matrix (c j7 ) is Ci, • • •, c m . Then by (18) each diagonal term 
is ± Ci, whence m =t= 3. In the rth and sth rows, ± c st lies above c sp , 

* Comptes Rendus, Paris, 104, 1887, 566-7. 

t TWorie des nombres, 1891, 294. 

t Assoc, frang. av. sc, 1896. 

§ L'interm^diaire des math., 3, 1896, 259-262. 

|| Periodico di Mat., 14, 1899, 125-139; also as a pamphlet of 1896. 
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and =F c sp above c st . Hence m must be even. It is assumed that, if 
a, = d = {% > m/2), (17) reduces to a Y^m square formula. Thus if 
m = 2*&>, where &> is odd, we would get an &> square formula by continued 
halving. Hence « = 1, m = 2 k . The impossibility of a 16 square for- 
mula is established more simply than in the earlier papers. 

27. K. Th. Vahlen* noted the error in Studnicka's deduction of the 
8 square theorem from the product of two two-rowed determinants and 
deduced that theorem by use of the product of two three-dimensional 
determinants [as had Antomari, §22] : 

(aa' + W + cc' + dd')(xx' + yy' + zz' + tt') 

= (ax + by + cz + dt)(a'x' + b'y' + c'z' + d't') 
+ (- b'x + a'y + dz' - ct'){- bx' + ay' + d'z - c't) 
+ (- c'x - dy' + a'z + bt')(- ex' - d'y + az' + b't) 
+ (- d'x + cy' - bz' + a't)(- dx' + c'y - b'z + at'). 

For a = a', etc., this gives the formula (4) for 4 squares. If a' is the con- 
jugate to a, it gives the 8 square theorem. He gave an analogous, much 
longer, formula which for a = a', etc., becomes the 8 square formula, but 
when a' is the conjugate to a, etc., does not yield a 16 square formula. 

28. E. Barbettef discussed the 4 and 8 square theorems in connection 
with magic squares. 

* Giornale di Mat., 39, 1901, 181^. 

t Les sommes de p-iemes puissances . . . , Liege, 1910. 



